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Abstract 



Using the NRQCD factorization formalism, we calculate the total cross 
section for the photoproduction of h c mesons. We include color-octet and 
color-singlet mechanisms as well as next-to-leading order perturbative QCD 
corrections. The theoretical prediction depends on two nonperturbative ma- 
trix elements that are not well determined from existing data on charmonium 
production. For reasonable values of these matrix elements, the cross section 
is large enough that the h c may be observable at the E831 experiment and at 
the HERA experiments. 
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Of all the charm-anticharm quark boundstates lying below the threshold of open charm 
production, the most elusive to experimental investigation is the h c meson. To date this 
particle has been observed by only a few experiments, which measure an average mass of 
M hc = 3.52614 ± 0.00024 GeV [Q. Studies of the h c are difficult because it has quantum 
numbers J PC = and thus cannot be produced resonantly in e + e~ annihilation, or appear 
in the decay of a J PC = 1 charmonium state via an electric dipole transition. 

Since charm quarks are heavy compared to Aqcd, it is natural to view charmonium as 
a nonrelativistic system, where the h c is a spin-singlet P-wave state of a c and c. This 
approach is taken in the color-singlet model |2[] (CSM) of quarkonium production and decay. 
In the CSM, it is assumed that the cc must be produced in a color-singlet state with the 
same angular-momentum quantum numbers as the charmonium meson which is eventually 
observed. However, the CSM has serious deficiencies. It is well known that perturbative 
QCD calculations of production and decay of P-wave quarkonia within the CSM are plagued 
by infrared divergences ||. More recently, measurements made at the Fermilab Collider 
Detector Facility (CDF) show that the CSM also fails to accurately predict the production 
cross sections of S-wave quarkonia (J/i/j and if)') [|J. 

The naive CSM has been supplanted by the NRQCD factorization formalism of Bodwin, 
Braaten, and Lepage |J. This formalism allows the infrared safe calculation of inclusive 
charmonium production and decay rates. It also predicts new mechanisms in which a cc 
pair is produced at short distances in a color-octet state, and hadronizes into a final state 
charmonium nonperturbatively. These color-octet mechanisms can naturally account for the 
CDF data on J /if) and if)' production. 

In this paper, we examine photoproduction of h c within the NRQCD factorization for- 
malism. Color-octet mechanisms play an essential role. By itself the CSM yields infrared 
divergent expressions for the h c cross section; however, once color-octet contributions are 
included we obtain sensible predictions. The resulting expression then depends on two un- 
determined nonperturbative matrix elements. Using heavy-quark spin symmetry these can, 
in principle, be estimated from similar matrix elements extracted from data on the produc- 
tion and decay of Xcj mesons. Currently, theoretical and experimental uncertainties in the 
determination of the Xcj matrix elements preclude us from making definitive predictions 
for the h c photoproduction cross section. This will be discussed in greater detail below. 
However, for some choices of the matrix elements which are consistent with current data, we 
find a large enough cross section that the h c may be observed in current photoproduction 
experiments at the DESY HERA collider and at the Fermilab fixed target experiment E831. 

In the NRQCD factorization formalism, inclusive quarkonium production cross sections 
have the form of a sum of products of short- distance coefficients and NRQCD matrix ele- 
ments. The short-distance coefficients are associated with the production of a heavy quark- 
antiquark pair in specific color and angular-momentum states. They can be calculated using 
ordinary perturbative techniques, and are thus an expansion in the strong coupling constant 
a s . The NRQCD matrix elements parameterize the hadronization of the quark- ant iquark 
pair, and each scales as a power of the average relative velocity v of the heavy quark and 
antiquark as determined by the NRQCD velocity-scaling rules 0. 

The NRQCD factorization formula for h c photoproduction, at leading order in v, is 

a ( 1 + N ^ hc + X) = JdxJ2fi/N{x) [a( 1 + i^cc(8, 1 S )+X)(O h 8 cCS ))^ + 
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a( 1 + i^ C c(l, l P x ) + X; riiOtCPi)) 



(1) 



where /i/jv(x) is the probability of finding a parton i in the nucleon with a fraction a; of the 
nucleon momentum. a(j + i — > cc(8, 1 5'o) is the short- distance coefficient for producing 
a cc pair in a color-octet 1 Sq configuration, and a(j + i —>■ cc(l, x Pi) + X; fi) is the short- 
distance coefficient for producing a cc pair in a color-singlet l P\ configuration. The matrix 
element (O x ^( 2S+1 Lj)) describes the hadronization of a color-singlet (color-octet) 2S+1 Lj 
cc pair into an h c . According to the w-scaling rules, both matrix elements in Eq. ([I]) scale as 
v 5 . The scale \x appearing in Eq. (|]) arises from the factorization of the cross section into 
long-distance and short-distance contributions. As we discuss below, the \i dependence of 
the color-octet matrix element is cancelled by the \x dependence of the color-singlet short- 
distance coefficient so that the expression for the cross section is \x independent. Note that 
in addition to the NRQCD scale p shown explicity in Eq. ([!]) there is a renormalization 
scale and a factorization scale (associated with the parton distribution function /i/jv(x)) 
which have been suppressed. These three scales do not have to be the same, however, in our 
numerical calculations we choose them to be equal. 

The short-distance coefficients in Eq. (^ can be calculated using the techniques of 
Ref. [0]. At leading order in a s , we obtain the following expression for the color-octet 
short-distance coefficient: 



da 



cc 



.8,%)) 
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4m5 



5(1 -z) 



(2) 



where z = p/x, p = im^/S^N, and S^n is the photon-nucleon center-of-mass energy squared. 
We calculate the color-singlet cross section, regulating the infrared divergences using dimen- 
sional regularization. The result in the MS scheme is: 



^( 1 + g^cc(l, 1 P 1 )+g;p) 



%n 2 aa 2 e 2 
27ml 
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where 



z 5 \n(z) z 2 (5 + 3z + 14^ 2 + 2z 3 + 9z 4 - z 5 ) \n(z) 



1 + z) 
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and the functional distribution is defined by 
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dz f(z) 
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The NRQCD expression for the cross section is obtained by substituting Eqs. (§) and @ 
into Eq. (|I|). Note that a 1/e divergence in the color-singlet coefficient has been absorbed into 
the definition of the leading color-octet matrix element. The renormalized matrix element 
(O^^Sq))^ can be shown to obey the renormalization group equation ||, 
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(6) 



Thus the logarithmic dependence on /i of the short- distance coefficient a (7 + g — > cc(l, 1 P\ + 
X; jj) is canceled to this order by the \x dependence of the renormalized matrix element 

Note that the leading color-octet coefficient is 0(a s ) while the leading color-singlet co- 
efficient is 0(aq). Therefore, next-to-leading order QCD corrections to the color-octet coef- 
ficient are of the same order as the leading color-singlet coefficient, and must be included if 
we are to have a complete 0(a 2 ) calculation of h c photoproduction. The 0(a 2 ) corrections 
to the color-octet contribution are computed in Ref. ||. These corrections are included in 
our calculation of the cross section. 

In order to make a prediction for the h c photoproduction cross section, we must deter- 
mine the values of the nonperturbative matrix elements (O^^Sq))^ and (O^^Pi)). To 
eliminate large logarithms in the expression for the cross section we choose /i = 2m c . At 
this time there does not exist a direct measurement of these matrix elements. However, they 
are related to similar matrix elements for Xcj production and decay by the (approximate) 
heavy quark spin symmetry of NRQCD: 

(OreS 1 )) = (O h s c( 1 S ))(l + O(v 2 )). (7) 

The size of the 0(v 2 ) corrections in Eq (0) can be estimated by studying radiative Xc and 
ip' decays, where predictions based on heavy quark spin symmetry agree with experiment 
to 20% accuracy ||. 

(Oi cl ( 3 Pi)) can be extracted from Xcj decays, or from the decay B — > Xcj + X. The 



authors of Ref. [10] calculate inclusive hadronic Xcj decay including next to leading order 



a s corrections. The result of their fit is: 

(OffPi)) 



m 2 



0.115 ± 0.016 GeV 3 . (8) 



The error includes only experimental uncertainties. It is particularly important to note that 
Ref. |10| does not include 0(v 2 ) relativistic corrections which are numerically of the same 



size as the 0(a s ) perturbative corrections that have been included. Relativistic corrections 



to the J /if) decay rate are large JTTJ. In the case of Xcj decay, relativistic corrections need 
to be analyzed before (Oi cl ( 3 Pi)) can be extracted with confidence. 

The decay B — > Xcj + X is calculated in Ref. [|I^] . Measurements of B decay []T3 1 allow 
an extraction of (Of cl ( 3 Pi)): 

(0fl( 9 3Pl)) = 0.42 ± 0.16 GeV 3 . (9) 

Again, the error quoted above is only due to experimental uncertainties. The authors of 
Ref. [0] state that their calculation suffers from large theoretical uncertainties due to next 



to leading order QCD corrections to the Wilson coefficients and the subprocess b — > ccs. 
Note that the two extractions of (0± cl ( 3 Pi)) agree only at the 2a level. 
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{Of 1 ( 3 5i)) can be extracted from the decay B — > Xcj + A, an d from CDF data on Xcj 
production. The result of the fit to B decay, after running from the scale m& to 2m c and 
converting from a cutoff regularization scheme to dimensional regularization, is 

{Of 1 ( 3 5 1 )) (2mc) = (2.9 ± 2.0) x 1(T 2 GeV 3 . (10) 

As before only experimental errors are included. The result of two different fits to Tevatron 
data are 

{Of^Sx)) = 0.98 ±0.13 x 10~ 2 GeV 3 (11) 
{Of 1 ( 3 ^i)) (2mc) = 2.6 x 10" 2 GeV 3 , (12) 



where the values are taken from Refs. [14] and [15| respectively. The error of Ref. [14 



includes only experimental uncertainty; Ref. [IS] does not quote errors. The central values 
of both fits lie within the la error of the extraction from P-decays. 

Note that in the Xcj production calculations of Ref. [|T3|| , there is no color-singlet con- 



tribution that gives an infrared divergence which needs to be absorbed into the definition 
of the color-octet matrix element, as is done in our calculation. Therefore, it is not pos- 
sible to relate the the "bare" matrix element appearing in Eq. (|TT|) with the renormalized 
matrix element needed for our calculation. However, the fragmentation calculation of Xcj 



production carried out in Ref. |15j makes use of the g — > Xcj fragmentation function. This 
fragmentation function includes both color-octet and color-singlet contributions, and has an 
infrared divergence in that is absorbed into the definition of the color-octet matrix element. 
This allows us to import the extracted value given in Eq. ([12]) into our calculation. 

Not all values for (0£ c (%)) and (O^Pi)) result in a physically sensible prediction for 
the cross section. Once the infrared divergence from the color-singlet contribution to the 
cross section is factorized, the remaining finite contribution is actually negative (for positive 
(Oj ic ( 1 Pi))/m 2 ). If the ratio of the color-octet to color-singlet matrix elements is too small, 
it is possible to obtain physically meaningless results. This is the case if, for example, the 
central values of the color-singlet matrix element extracted from P-decays (Eqs. (0,110)) is 
used. Therefore it is impossible to put a lower bound on the cross section given our current 
state of ignorance concerning this matrix element. It is also important to point out that 
the NRQCD velocity scaling rules imply that (0£ o ( 1 Sb)) (am ° ) and {O^ i 1 Pi)) / m 2 c should be 
roughly the same size. Thus, theory would prefer a smaller value for the matrix element 
(Of c ( 1 P 1 )), as suggested by the fit to Xcj decay. Clearly more accurate extractions from 
the Tevatron, from P-decays, and from Xcj decays are needed to clarify the situation. This 
may be possible once next-to-leading perturbative QCD corrections and leading relativistic 
corrections to these processes are calculated. 

For the NRQCD matrix elements, we will use {0\ c : ( 1 Pi)) / 'm 2 = 0.115 GeV 3 as de- 
termined by the analysis of Xcj decays (Eq. (§)), and we choose (0g c ( 1 ,So)) < ' 2mc - ) = 2.6 x 
10 _2 GeV 3 . We use the CTEQ 3M parton distribution function with the factorization scale 
choosen to be equal to the NRQCD scale fi = 2m c . The resulting cross section is plotted 
as the solid line in Fig. |I[ However, our results are extremely sensitive to the uncertainty 
in the determination of the matrix elements. To show this, we also plot, as the dashed line, 
the cross section with the choice (O^^Pi)) /m 2 = 0.2 GeV 3 . The cross section then drops 
by roughly a factor of four. 
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FIG. 1. The total h c photoproduction cross section 17(7 + N — > h c + X) as a function of the 
center of mass energy ^/S 7 n- The solid curve is for the choice (O^ ( l Pi)) / m 2 c = 0.115 GeV 3 and 
(Og c ( 1 5o)) (2mc) = 2.6 x 10~ 2 GeV 3 . The dashed line is obtained by keeping the same value for 
the color-octet matrix element and by changing the value of the color-singlet matrix element to 
(O^^P^/ml = 0.2 GeV 3 . 
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We conclude with a brief discussion of the possibility of observing the h c . The E831 
photoproduction experiment corresponds to roughly J S^n — 20 GeV, while the HERA 



experiments take photoproduction data at approximately J S^n — 100 GeV. At these 
energies the cross section for h c production, assuming our first choice of NRQCD matrix 
elements, is 30 nb and 62 nb respectively. These cross sections are comparable to J/ip 
production. The h c can be detected either via the rare decay h c — > J /ip + 7r with branching 
ratio on the order of 1% [|T(J, or the radiative decay h c — > r] c + 7 (BR = 50%) |T7|] . If it is 



possible to reconstruct these decay modes, then the possibility of observing the h c is real, 
and should be experimentally investigated. 
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